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Abstract —Nuclear energy plays a key role in long-term development plans and can guarantee the supply of electricity to 
some regions. On the other hand, the implementation of these projects tends to require long maturation periods, require high 
investment costs and may be a source of pollutants, such as Lead (Pb). For this reason, the periodic environmental 
monitoring of the concentration of pollutants becomes necessary, according to the current legislation. Thus, the present work 
presents as a proposal a method of optimization of Pb quantification in environmental analyzes. The study was developed 
through the monitoring of wastewater samples from INB - Industrias Nucleares do Brasil. From the Fractional Factorial 
Design(FFD) the most representative variables of the Pb isolation process were determined, and through the Central 
Composite Design (CCD) the response surface was found, generating a regression model that represents the system. The 
results indicated that the optimization of the chemical yield of Pb is associated to the higher dosages of Nitrilotriacetic Acid 
(NTA) and Sulfuric Acid (H2SO4) in the samples. The optimum yield condition was obtained in the region of200 ml ofH2S04 
and 4.0 g of NTA, considering the range tested. The proposal proved to be effective for the validation tests of the model, 
obtaining an increase of up to 32% in the Pb yields of the analyzes. 

Keywords—Chemical Yield, Environmental Analyzes, Lead, Nuclear Power, Response Surface. 

I. Introduction 

The generation of electricity from nuclear sources can be considered advantageous and represent a viable strategy of 
diversification of the energy matrix in some countries [1,2], meeting the goals of reducing greenhouse gas emissions 
compatible with the Kyoto Protocol [3,4]. Nuclear energy plays a key role in long-term development plans and can guarantee 
energy supply for some regions [5]. On the other hand, the implementation of these projects tends to require long maturation 
periods, from the installation of the plant until its operation, in addition to requiring high investment costs and representing 
sources of environmental degradation [4,6,7]. 

Currently in Brazil nuclear sources have little expressiveness in the electric matrix, but this scenario tends to undergo 
changes until the year 2050 [8]. The Brazilian government plans to build new nuclear plants by 2030 [6,9], which should 
represent an increase of approximately 6 GW in capacity country [6,10]. 

The process of generating energy from nuclear sources occurs through nuclear fission in reactors, which are fed by fuels from 
the extraction of uranium ore [11,12,13,14,15]. From the mining stage to the production of nuclear fuel, innumerable 
pollutants from soil, air and water, especially heavy metals, such as lead, which has high toxicity and may cause serious 
social and environmental impacts [16,17,18,19,20,21,22], 

In the case of wastewater from nuclear activities, there is a demand for specific treatments so that they can meet the quality 
adjustments [23,24], as established by Resolution No. 430 of the Conselho Nacional de Meio Ambiente - CONAMA [25], 
According to this resolution, release of effluents containing Pb to water bodies or collection systems must comply with the 
maximum acceptable limit of 0,5 mg.F' 1 [25]. Therefore, periodic environmental analyzes must be carried out to monitor and 
control the concentration of Pb according to current legislation [26,27,28,29]. Such analyzes involve a methodology of 
isolation and quantification, having as principle the separation of the lead from the other elements, besides the separation 
between them by selective precipitation aiming the determination of the yield of Pb in the samples. 

Thus, the present study proposes the optimization of the lead insulation process in environmental analyzes, as a way to 
improve the efficiency of the quantification of Pb yields, in order to guarantee reliable results that comply with the norms of 
the current legislation. The Fractional Factorial Design (FFD) was used to determine the most representative variables of the 
process [30,31,32,33,34], in order to verify the relevant sources of interference of the method employed. The Central 
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Composite Design (CCD) was used to find the response surface where the yield of Pb is maximum, generating a regression 
model that represents the system [35,36,37,38,39,40]. 

II. Material and Methods 

The execution of the factorial design was started with the collection of a sample composed of water (60 liters), collected from 
several monitoring points of the company Industrias Nucleares do Brasil (INB), a unit of Caldas-MG. The sample was 
filtered and reduced in hotplate up to 30 liters. In the first step of the optimization, four variables were selected to be tested in 
the chemical insulation of lead (Pb): i) addition order of Nitrilotriacetic Acid (NTA), Sodium Hydroxide (NaOH) 6 mol.L' 1 , 
Ammonium Sulfate (NH 4 SO 4 ) 25 mg.L " 1 and heated distilled water; ii) NTA dosage; iii) dosage of Sulfuric Acid (H 2 S0 4 ) 3 
mol.L 1 , and iv) dosage of Sodium Sulphide (Na 2 S) 1 mol.L' 1 . 

According to the Fractional Factorial Design procedures, the number of experiments to be performed is of the order of 2 n l , 
where n represents the number of variables listed [30]. Thus, in this study, 8 experiments (2 41 ) were performed in triplicate, 
totaling 24 tests. The FFD codifies the possibilities to be tested in the level tests. In this case, the level -1 represents the 
currently applied conditions and the level +1 the new conditions tested. Table 1 shows the cataloged variables and their 
respective levels. 


Table 1 

Variables listed and their respective levels. 


Variable 

Level 

A) Order of addition of the reagents: nitrilotriacetic 
acid (NTA), sodium hydroxide (NaOH) 6 mol.L" 1 , 
ammonium sulfate (NH 4 S0 4 ) 25 mg.L " 1 and 
heated distilled water 

Level (-1): 1° NTA, 2° NaOH, 3° NH 4 S0 4 , 4° heated distilled 
water (current order). 

Level (+1): 1° heated distilled water, 2° NaOH, 3° NTA, 4° 
NH 4 S0 4 (new order). 

B) Dosing of NTA 

Level (-1): 2,0 g of NTA (current dosage). 

Level (+1): 4,0 g of NTA (new dosage). 

C) Dosing of H 2 S0 4 3 mol.L ' 1 

Level (-1): 50 ml of H 2 S0 4 0,1 mol.L " 1 (current volume). 

Level (+1): 150 ml of H 2 S0 4 0,1 mol.L ' 1 (new volume). 

D) Dosing of Na 2 S 1 mol.L " 1 

Level (-1): 1,0 ml of Na 2 S (current volume). 

Level (+1): 3,0 ml of Na 2 S (new volume). 


The next step consisted in inserting all variables and their respective levels into Minitab software (version 16.2.4.4), in order 
to obtain an array with all possible combinations of experimental tests (Table 2) [41], 

Table 2 


Experimental matrix of the Fractional Factorial Design for ] 

Pb. 

Experiment 

Level 

NTA (g) 

H 2 S0 4 (ml) 

Na 2 S(ml) 

1 

+ 1 

2 

50 

3 

2 

-1 

2 

50 

1 

3 

-1 

4 

150 

1 

4 

+ 1 

4 

150 

3 

5 

+ 1 

4 

50 

1 

6 

-1 

2 

150 

3 

7 

+ 1 

2 

150 

1 

8 

-1 

4 

50 

3 


From the matrix obtained, the tests of isolation and quantification of Pb with the 24 samples were performed, following the 
analysis procedures described below; 
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1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 


9. 

10 . 

11 . 

12 . 


13. 


14. 


15. 


After the filtration initially described, the sample was transferred to 24 beakers of 2 liters, each containing 1 liter of 
sample; 

Then 1 mL of 20 mg.mL " 1 Barium solution and 1 mL of 20 mg.mL " 1 Lead solution were added with a micropipettor; 

With dispensing, 5 mL of Citric Acid solution (C 6 H 8 0 7 ) 1 mol.L " 1 and 2 drops of 0.1% Methyl Red in each beaker were 
added. The pH of the sample was adjusted between 6.0 and 7.0 with Ammonium Hydroxide (NH 4 OH) and Nitric Acid 
(HN0 3 ) 5 mol.L ' 1 slowly, until turning from red to yellow, using calibrated Phmetro; 

To each beaker were added 50 ml of Sulfuric Acid solution (H 2 SO 4 ) 3 mol.L" 1 , with stirring, and the solution was 
heated to boiling and allowed to settle until the following day; 

The supernatant was then drained by siphoning and the precipitate transferred quantitatively into the glass centrifuge 
tube using osmosis purified water; 

The sample was centrifuged for 10 minutes at 2000 RPM and the supernatant discarded by siphoning. This was done 
with 20 ml of H 2 S0 4 solution 0.1 mol.L " 1 added to the precipitate, while the supernatant was centrifuged and discarded 
by siphoning again; 

To the precipitated was added 2.0 g of Nitrilotriacetic Acid (NTA), 10 mL of Sodium Hydroxide solution (NaOH) 6 
mol.L" 1 , 50 mL of distilled water and 5 mL of Ammonium Sulfate (NH 4 ) 2 S0 4 ) 25 mg.mL" 1 ; 

The tubing solution was heated in a water bath at approximately 90 °C for 15 minutes until the precipitate dissolved 
completely. Then, the pH of the sample was adjusted by adding glacial Acetic Acid to pH between 4.0 and 4.5 and the 
solution was held until the following day; 

The solution obtained was centrifuged, the supernatant was transferred to a 250 mL beaker to separate it from the 
precipitate and it was reserved in tubes; 

In each of the supernatants obtained was added 1 mL of Sodium Sulfate (Na 2 S) solution 1 mol.L" 1 . These were then left 
in a water bath at 90 °C for approximately 15 minutes; 

After cooling the samples in the chambers, the solutions obtained were centrifuged, the supernatant was discarded by 
siphoning and the precipitate was washed with 50 ml of distilled water. The sample was then centrifuged again and the 
supernatant discarded by siphoning. To the precipitate was added 10 mL of HNO 3 solution 5 mol.L" 1 ; 

The solution was taken in a water bath at approximately 90 °C and waited from 100 to 120 minutes until the total Lead 
Sulphide formation reaction occurred, which is indicated by the absence of dark, sulfur-forming spots. Each of the 
samples was screened on white strip paper and the filter was washed 3 times with heated distilled water. They were then 
left in a water bath for 10 minutes; 


The pH of the resulting solution was adjusted between 4.0 and 4.5 with Ammonium Acetate solution (NH 4 CH 3 COO) 
40%. To each sample was added 1 mL of Bismuth solution 5 mg.mL " 1 and 2 mL of Sodium Chromate solution 
(Na 2 Cr0 4 ) 30%; 

The material was filtered using a vacuum filtration kit, with 0.45 pm membrane filter and diameter 0 = 47 mm, 
previously weighed in analytical balance. The precipitate was washed with distilled water, the filter was removed from 
the Kit and stoved at 60 °C for 90 minutes; 

Finally, the filter was removed from the oven and placed in desiccator for 10 minutes. The samples were taken from the 
desiccator, weighed and the yield was calculated considering the added standard mass equal to 31 mg. The gravimetric 
yield was calculated from the Equation: 


Y grav (%) = 100^ 


(1) 


Where: M p is the mass of the filter with precipitate; 

M f is the mass of the filter; 

M a is the theoretical mass of Lead standard (31 mg). 


16. After the previous tests were carried out to determine the response surface of the CCD, the filters were digested in a 
microwave oven and then analyzed in an Inductively Coupled Plasma Atomic Emission Spectrometer (ICP). The 
chemical yield was calculated from the Equation: 


X:hem C^ 0 


c Pb -Vf 


1 Pb ■ m PbCrO 4 


.100 


( 2 ) 


PbCrO 4 
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Where: C Pb is the concentration of Pb (g.L 1 ); 

V f is the volumetric flask volume (0.05 L); 

MMp b is the molar mass of Pb (207.19 g.mol' 1 ); 
nipbcro 4 is the mass of Lead Chromate (0.031 g); 

MMpbCro 4 is the molar mass of Lead Chromate (323.2 g.mol’ 1 ). 

2.1 Validation of the model 

To validate the model found by the response surface of the Central Composite Design (CCD), new tests were performed for 
four different samples of wastewater under the current conditions (50 mL of PLSCL and 2 g of NTA) and under the optimum 
conditions established by the model. The assays followed the same procedure described previously for the chemical isolation 
of Pb by selective precipitation. 

III. Results and discussion 

After the Pb isolation for the 24 samples of FFD, the average gravimetric yields were obtained in each of the experimental 
conditions (Equation 1). In the next step, the variables were coded at levels -1 and +1, as described in Table 1, obtaining a 
new coefficient matrix. From the possible combinations between the four variables selected in the coded levels, the contrasts 
involved in the design and their respective percentages of contribution were calculated (Table 3). 

Table 3 


Contrasts calculated for the development of the Fractional Factorial Design for Pb. 


Combinations 

A+BCD 

B+ACD 

C+ABD 

D+ABC 

AB+CD 

AC+BD 

AD+BC 

Contrasts 

-5,88 

-3,88 

1,04 

-2,29 

1,79 

0,04 

-1,29 

Contrasts (%) 

56,82 

24,72 

1,78 

8,64 

5,28 

0,003 

2,746 


To verify which of the selected variables have a significant effect on the system, the confidence interval of the contrast error 
was calculated. First, the number of degrees of freedom (DF) and the quadratic sum of the pure error (QSPE) were 
calculated. From the DF and QSPE values we found the experimental variance (S 2 ), the experimental error (S), the contrast 
error (S c ) and finally the confidence interval of the contrast error (I Sc ) for a 95% confidence interval, as shown in Equations 
(3), (4), (5) and (6), respectively. All statistical data obtained are presented in Table 4. 


c 2 _ £ GL . SQEP 

~ Z GL 

(3) 

s = Vs 2 

(4) 

s c = — 

c 3 

(5) 

Is c =t.S c 

(6) 


Table 4 

Statistical data found for Pb analysis 


Experimental Variance (S 2 ) 

10,20 

Experimental Error (S) 

3,20 

Contrast Error (S c ) 

1,80 

Confidence Interval of the Contrast Error (I Sc ) 

3,30 


As the confidence interval for the Contrast Error found was 3.30, only the contrasts larger than this value in modulus will be 
significant. Thus, as shown in Table 4, only the interactions A+BCD (contrast: -5.88) and B+ACD (contrast: -3.88) are 
relevant. 
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As the probability of occurrence of a third-order interaction is minimal, it is considered that only A and B are influencing the 
analysis, only the dosage of H 2 SO 4 (56.8%) and NTA (24.7%). According to researches by Ozer and Gu£er (2011) [42], the 
volume of Sulfuric Acid presented as a relevant variable in the digestion process of Pb samples, since it can participate in the 
formation of Lead Sulphate precipitate. Studies by Sakurai et al. (2006) [43] also reported the significant effect of H 2 S0 4 on 
the plastic digestion for the determination of Pb in a microwave oven. No results were found for NTA in the literature. 


To visualize what was found by the previous calculations, a graph of cubes was constructed (Fig. 1), with the aid of the 
software OriginPro 9.1.0, which presents the variables under study and their respective levels coded as a function of the 
obtained Pb yield [44]. 


Cube Plot (data means) for Rendimento Pb 




-1 1 

H2S04 

-1 1 
Na2S 


Figure 1. Cube plot for Pb yield. 


Fig. 1 shows that level -1 (current dosage) for FLSO 4 has a higher average yield than level +1 (new dosage) for the same 
reagent, thus contributing to the maximization of Pb yield at level -1. For the NTA dosage, the level -1 also has an average 
yield greater than the level +1, evidencing that the level -1 is more advantageous for the increase of the Pb yield. 

3.1 Response surface 

The initial proposal to define the response surface of the FFD in question was to employ a Central Composite Design (CCD) 
with two variables (dosage of Sulfuric Acid and NTA). This type of experiment employs compounds with a central point in 
order to identify the best working conditions of an analytical method. The variables are tested in 5 different levels, ranging 
from W2 to a/2 [30], 

By the established method, the coded interval (Ic) should be 2.83 for both H2SO4 and NTA dosage. The actual interval (Ir) 
was defined as 200 ml for H 2 S 0 4 and 3.5 g for NTA. From the values of Ic, Ir and the levels coded according to the Central 
Composite Design were calculated the actual levels. Table 5 shows the actual and coded levels defined by the CCD 
technique. 
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Table 5 


Coded and actual levels established by the Central Composite Design method. 


Experiment 

Coded Levels 

Real Levels 

Sulfuric Acid 

NTA 

Sulfuric Acid 

NTA 


3 mol.L' 1 

3 mol.L' 1 (ml) 

(g) 

1 

-1 

-1 

29,3 

1,0 

2 

0 

0 

100,0 

2,3 

3 

0 

-1,4142 

100,0 

0,5 

4 

-1 

1 

29,3 

3,5 

5 

-1,4142 

0 

0,0 

2,3 

6 

1 

-1 

170,7 

1,0 

7 

0 

1,4142 

100,0 

4,0 

8 

1 

1 

170,7 

3,5 

9 

1,4142 

0 

200,0 

2,3 


After the experiments were carried out by the method described for the lead insulation, the filters were obtained in a 
microwave oven and afterwards an Atomic Emission Spectrometer with Inductively Coupled Plasma (ICP). This procedure 
was performed in order to verify if the filters only have Lead Chromate (resulting from the analysis), or if any interfering 
substances are still present that were not totally removed in the isolation procedure and are contributing to the different 
gravimetric yield got from the chemical, obtained by quantification of Pb in the ICP, masking the true result of the analysis. 

After the tests, the gravimetric yields were calculated by the Equation (1) and the chemical yields by Equation (2) of the 
methodology. The results indicated differences ranging from 5.6 to 33.6% among the yields obtained by gravimetry and 
quantification of Lead. The systematic difference observed, always with the gravimetric yield greater than the chemical, 
made it possible to list some possibilities on the origin of the difference, such as the addition of Bismuth solution 5 mg.L" 1 in 
the Pb insulation and the residual moisture contained in the due to their drying temperature (60 °C), which may be 
contributing to a larger mass in the filter. Because the results obtained by quantifying Pb in ICP are more accurate than those 
found only by gravimetry, it was decided to use the results of chemical yield for the realization of the CCD. 

With the execution of the method, it was verified that only the variables NTA and the interaction Sulfuric Acid*NTA are 
significant for the model. In this way, the other variables and interactions were extracted, being possible to obtain in the 
software Minitab 16.2.4.4 the following response surface model, which results in the conditions that provide the highest 
yields [41], 

Pb Yield (%) = 63,6828 + 3,7383. X 1 + 3,4385. (X^) (7) 

Where X) and X 2 are the coded levels of NTA and Sulfuric Acid, respectively. 

Analyzing the model obtained, the greatest influence of the NTA variable (higher coefficient) is verified, although the 
interaction between the variables was also significant, evidencing that the response presents a distinct tendency, according to 
the established levels. Both variables show a positive coefficient, indicating that the higher the levels, the higher the yield. It 
is noteworthy that Sulfuric Acid alone was not significant. 

Studies by Ozer and Gu 5 er (2011), Savic et al. (2014) and Zare-Dorabei et al. (2016), involving Pb analyzes, also 
demonstrated the direct proportionality relation in their response surface models, since the variables were significant and 
positive coefficients, therefore, the higher their levels, the greater the associated response [42,45,46]. The response surface 
graph for the maximization of the chemical yield of Pb was elaborated as a function of the NTA and H 2 S0 4 dosage, through 
the OriginPro 9.1.0 software, as shown in Fig. 2 [44]. 
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51,50 

53,94 

56,38 

58,82 

61,26 

63,70 

66,14 

68,58 

71,02 

73,46 

75,90 


Figure 2. Response surface obtained for the chemical Lead yield as a function of the NTA 

and Sulfuric Acid 3 mol.L" 1 dosages. 

Fig. 2 shows that higher concentrations of H 2 SO 4 and NTA lead to higher yields, so the optimum condition in which the 
highest chemical yield of Pb for the range tested would be obtained in the region of 200 ml of H 2 SO 4 and 4.0 g of NTA. It is 
worth mentioning that the current methodology makes the dosage of 2 g of NTA and 50 mL of Sulfuric Acid, which 
corresponds to a yield of about 62%. 

3.2 Validation of the model 

The results found for the validation tests of the response surface model for the current dosage (50 ml of H 2 S0 4 and 2 g of 
NTA) and for the new dosage (200 ml of H 2 S0 4 and 4 g of NTA) showed that after applying the optimum conditions 
established by the response surface model obtained there was an increase of up to 32% in the chemical yields of those 
samples. These results can validate the model found and confirm the significant improvement that the application of the 
Factorial Design provides to the chemical isolation of Lead and optimizes the yields resulting from the analyzes. 

IV. Conclusion 

> The results of the optimization of the chemical isolation of Pb through the FFD showed that only the H 2 S0 4 dosage and 
the NTA dosage of the isolation methodology influenced the analysis; 

> By performing the CCD to obtain the response surface, it was verified that only the variation in the NTA dosage and in 
the NT A/Sulfuric Acid interaction were significant; 

> After the exclusion of the non-significant variables, the response surface showed that the higher the NTA and H 2 S0 4 
dosage, the higher the associated chemical Pb yield. The optimal condition at which the highest chemical yield of Pb is 
obtained for the range tested is the region close to 200 ml of H 2 S0 4 and 4.0 g of NTA; 

> The variable H 2 S0 4 dosage alone was not significant in obtaining the response surface, only the NTA/H 2 S0 4 
association showed influence in the analysis; 

> Validation tests of the response surface model showed an increase of up to 32% in the chemical yields of the samples. 
These results confirm the optimization of the yields that the application of the Factorial Design provided to the chemical 
isolation of Pb. 
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